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1. Summary 

Flight-measured high Reynolds number turbulent- 
flow pressure distributions on a transport wing in 
transonic flow are compared to unstructured-grid 
calculations to assess the predictive ability of a three- 
dimensional Euler code (USM3D) coupled to an inter- 
acting boundary layer module. The two experimental 
pressure distributions selected for comparative analy- 
sis with the calculations are complex and turbulent but 
typical of an advanced technology laminar-flow wing. 
An advancing front method (VGRID) was used to 
generate several tetrahedral grids for each test case. 

Initial calculations left considerable room for 
improvement in accuracy. Studies were then made of 
experimental errors, transition location, viscous 
effects, nacelle flow modeling, number and placement 
of spanwise boundary layer stations, and grid resolu- 
tion. The most significant improvements in the accu- 
racy of the calculations were gained by improvement 
of the nacelle flow model and by refinement of the 
computational grid. 

Final calculations yield results in close agreement 
with the experiment. Indications are that further grid 
refinement would produce additional improvement but 
would require more computer memory than is 
available. 

The appendix data compare the experimental 
attachment line location with calculations for different 
grid sizes. Good agreement is obtained between the 
experimental and calculated attachment line locations. 

2. Introduction 

2.1. Background 

Code validation by comparison with flight test 
results is an important pent of determining the accu- 
racy and usefulness of a newly developed code. These 
comparisons can be very important in practical appli- 
cations such as determining manufacturing tolerances 
of turbulent-flow aircraft or designing an entirely new 
laminar flow aircraft configuration. Industry users 
urge code developers to deliver user-friendly products 
tested and validated for everyday use. 


Between 1986 and 1991, NASA and its industry 
partners completed a series of flight tests on various 
airplanes modified to obtain laminar flow (refs. 1 
through 7). These tests culminated in the successful 
flight validation of the hybrid laminar-flow control 
(HLFC) concept at transonic speeds and high 
Reynolds numbers on both wings and nacelles (refs. 2 
through 7). This concept involves a combination of 
active laminar-flow control (LFC) with suction from 
the leading edge to the front spar and passive laminar 
flow obtained by means of a favorable pressure gradi- 
ent from the front spar to more rearward locations. 
During these experiments, a high-fidelity database was 
created that is useful for the validation of advanced 
computational fluid dynamics codes. The flight data 
used in this report were obtained on a modem trans- 
port airplane with high-bypass-ratio wing-mounted 
engines. One result of the flight experiments was the 
realization that the design procedure used for the wing 
surface pressure computations was often inadequate 
for predicting details of a complex pressure distribu- 
tion. It was therefore decided to undertake a detailed 
study of the ability of an advanced aerodynamic code 
to predict three-dimensional experimental pressure 
distributions. 

2.2. Objective 

Flight-measured high Reynolds number turbulent- 
flow wing pressure distributions are compared to 
unstructured-grid calculations to assess the predictive 
ability of a three-dimensional Euler code (USM3D) 
coupled to an interacting boundary layer module. The 
transonic test points selected for analysis have pres- 
sure distributions that are complex and turbulent, yet 
typical of an advanced technology laminar-flow wing. 

2.3. Technical Approach 

The state of the art in computational fluid dynam- 
ics offers a variety of methods for computing aircraft 
aerodynamic performance. These methods range from 
simple linear methods to full Navier-Stokes simula- 
tions. When viscous effects are important (e.g., for a 
laminar- flow airplane), they must be included in the 
mathematical model used to compute the flow. In 
addition, the engineer must decide whether to use a 
structured or unstructured computational grid. The 
generation of grids suitable for computation of an 
inviscid flow is generally easier for unstructured-grid 



methods, while structured-grid methods are generally 
more efficient in terms of both the computer memory 
requirements and the execution time of the flow 
solver. 

The geometry in the present study is sufficiently 
complex to make generation of a block-structured grid 
costly and time-consuming. Any gains in solver effi- 
ciency would be small compared to the time required 
to generate a suitable grid. Such methods, however, 
are relatively mature and robust and efficiently pro- 
vide viscous solutions once the grid is generated. 

While unstructured-grid Navier-Stokes codes 
exist, the generation of unstructured grids suitable for 
high Reynolds number flows remains difficult. The 
grid must be highly stretched near solid boundaries to 
properly resolve the boundary layer without generat- 
ing an unnecessarily large number of grid points in the 
streamwise and crossflow directions. Methods have 
been developed to produce grids of this type, but at 
present they lack the robustness required for wide- 
spread use. 

The coupling of an inviscid unstructured-grid 
method with an interacting boundary layer model 
offers a good compromise for the present study. Invis- 
cid unstructured-grid methods and grid generation 
techniques are well established and relatively robust. 
When coupled to a reasonable boundary layer model, 
an inviscid solver will capture most features of interest 
in an attached flow. For these reasons, the present 
work uses this method to analyze the aerodynamics of 
the configuration. 

3. Nomenclature 


C L 

lift coefficient 

c p 

surface pressure coefficient 

H 

altitude 

HLFC 

hybrid laminar flow control 

IBL 

interactive boundary layer 

LE 

leading edge 

LFC 

laminar flow control 

M 

free-stream Mach number 


MAC mean aerodynamic chord 

R free-stream Reynolds number based on 

mean aerodynamic chord 

WBL wing butt line 

x/c chordwise coordinate, fraction of local 

chord 

a angle of attack, deg 

4. Overview of Experiment 

The following sections present a description of the 
overall experiment. One wing of a modem high- 
bypass-ratio transport airplane was modified by 
removing the two slats closest to and outboard of the 
nacelle. The slats were replaced with a new suction 
leading-edge test article and a Krueger flap. An exten- 
sive amount of experimental data was then obtained at 
many different flight test points that are representative 
of the operating conditions of a commercial airplane. 

4.1. Test Article Description 

Figure 1 shows a sketch of the perforated leading- 
edge test article with a sweep of approximately 25° 
that is mounted on the left wing of the airplane. The 
test article was 22.4 ft long and extended 25 in. from 
the leading edge to the front spar. Titanium sheet of 
0.040-in. thickness was used for the perforated outer 
surface of the suction panel. 

Four rows of 34 pressure taps each were installed 
outboard of the nacelle at wing butt lines (WBL’s) 
344, 389, 428, and 479 in. (fig. 1). Each row included 
20 upper surface taps mounted in the leading-edge box 
suction panel. The taps did not penetrate the perfo- 
rated titanium skin but terminated in a small cavity in 
the stringer beneath the skin. The resulting pressures 
measured on the leading-edge test article were thus 
relatively disturbance free. The 10 surface pressure 
measurements made in the sparbox region used strip- 
a-tube belts bonded to the wing surface. The forward 
edges of the belts were tapered to a wedge shape to 
reduce any flow disturbance. The most forward inspar 
pressure tap (e.g., 20 percent chord at WBL = 344 in.) 
was close to the initial wedge surface so that its pres- 
sure values might have been affected by flow accelera- 
tion over the wedge. A tapered plastic strip was 
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applied at the belt sides to help reduce flow distur- 
bances further. The pressure belts terminated at the 
rear spar. 

4.2. Test Cases 

Two transonic high Reynolds number test cases 
were chosen for data analysis. The first had a more 
forward shock location with M = 0.793 and 
a = 2.598° ( C L = 0.447). The second had a more 
aft shock location with M = 0.820 and a = 2.377° 
{C L = 0.449). Additional test characteristics for the 
two flight conditions are summarized in table 1. 

The suction system was not used when measuring 
the pressures over the wing box. Thus, the flight data 
represent a turbulent wing flow in both cases. 

5. Overview of Calculations 

More detailed information concerning the compu- 
tational method used for analysis of the configuration 
is presented in the following sections. First, the fea- 
tures of the flow solver and boundary layer model are 
described. Next, the method used to generate grids is 
outlined. Finally, some issues pertaining to the air- 
plane geometry are discussed. 

5.1. Code Description 

The flow solver, referred to as a three-dimensional 
Euler code (USM3D) (refs. 8 through 10), computes 
the steady-state solution to the Euler equations on a 
grid of tetrahedral cells. Averages of conserved quan- 
tities (mass, momentum, and energy) are stored at the 
centroid of each cell. A cell-centered finite-volume 
formulation of Roe’s flux difference splitting scheme 
is used to discretize the spatial derivatives. A red- 
black Gauss-Seidel scheme is used to integrate the 
equations in time to the steady-state solution. 

To complete the specification of the solver, the 
boundary conditions used for the present calculations 
are given. On the wing, the normal velocity compo- 
nent is specified on the surface. On the fuselage, 
nacelle, pylon, and symmetry plane, flow tangency is 
enforced. In the far field, Riemann invariants are used 
to compute conserved quantities on the boundary. On 
the inlet and exit faces of the nacelle, conserved quan- 


tities are prescribed so that mass flow is conserved, 
and total pressure and total temperature within the 
simulated engine remain constant at free-stream con- 
ditions. The resulting nacelle flow is equivalent to that 
through an open duct. 

5.2. Boundary Layer Model 

The boundary layer on the wing is modeled as a 
set of two-dimensional strips at specified span loca- 
tions. At each of these boundary layer stations, the sur- 
face pressure computed by USM3D is used as the 
boundary condition in a two-dimensional boundary 
layer code (ref. 11). The solution of the boundary layer 
code yields the displacement thickness that is used to 
compute a transpiration velocity at the surface. The 
computed transpiration velocities are bilinearly inter- 
polated to the wing surface nodes between each pair of 
boundary layer stations and are then used as a bound- 
ary condition for USM3D. Input parameters allow 
specification of the number of flow solver iterations 
between boundary layer calculations. 

The boundary layer code uses the Keller box 
method to solve a finite-difference discretization of 
the two-dimensional compressible boundary layer 
equations (ref. 11). By using the velocity and density 
at the edge of the layer, the density at the wall, and the 
displacement thickness computed by the boundary 
layer code, an equivalent transpiration velocity is 
computed via Lighthill’s equation (ref. 12). 

5.3. Computational Grids 

The grid generation code used is referred to as 
VGRID (refs. 13 through 15) and is based on the 
advancing-front technique. The distribution of tetrahe- 
dra is governed by user-prescribed sources that influ- 
ence the solution of a Poisson equation on an 
underlying Cartesian background grid. The placement 
of sources on this background grid allows clustering of 
cells near points of interest. The process begins with a 
triangulation of each surface patch of the geometry. 
From this triangulation (which forms the initial 
“front”), points are added to form cells according to 
the spacing specified on the background grid. 

Several grids were generated during the course of 
this study, and the relevant parameters of each grid are 
summarized in table 2. Included in these parameters is 
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the grid spacing specified at various chord locations in 
the midspan region. The sequence of refinements rep- 
resented by grids 1, 3, and 4 in table 2 is illustrated in 
figure 2, in which each surface grid is depicted in the 
region where the experimental data were measured 
(i.e., the midspan region of the wing). The chord 
regions targeted for refinement were chosen based on 
the experimental measurements. Illustrations of the 
triangulation of the entire airplane surface for the fin- 
est grid (grid 4 in table 2) are presented in figure 3. 
The overall size (i.e., number of tetrahedra) of the fin- 
est grids was constrained to be less than about 1 .4 mil- 
lion cells, due to the memory limit (256 million 
words) of the computer when the implicit integration 
scheme is used. 

The implementation of the coupling between 
USM3D and the interactive boundary layer model 
allows boundary layer calculations only on chordwise 
patch boundaries on the wing. The initial grid (grid 1 
in table 2) had only 10 such boundaries along the span, 
as shown in figure 4, along with the spanwise loca- 
tions of the experimental data. Early results suggested 
that more stations may be required, so each wing patch 
was split in the chordwise direction, resulting in a 
maximum of 1 8 boundary layer stations, as shown in 
figure 5. 

5.4. Geometry 

A pointwise description of the aircraft surface was 
supplied by the manufacturer. Because the wing struc- 
ture was designed to match the specified geometry 
under a cruise loading, this geometry includes 
aeroelastic deflections of the wing. 

The initial aircraft geometry exhibited some sur- 
face irregularities on a few of its surface patches. 
These patches were located on the pylon and on the 
fuselage near the wing leading edge (LE). Efforts were 
made during grid generation to smooth out this effect. 
These regions are not expected to have a significant 
impact on the results in the region of interest. 

6. Results and Discussion 

Figure 6 presents a pressure distribution typical of 
a laminar flow control wing. The relative complexity 
of this pressure distribution makes it a good test case 
for an advanced three-dimensional code. The desired 


characteristics of the wing’s pressure distribution 
include an initial rapid flow acceleration, a slight peak 
in the pressure distribution near the leading edge, fol- 
lowed by a long, mild acceleration and then an aft 
pressure recovery region. For applications such as cor- 
rectly predicting the amount of laminar flow obtained 
over a wing’s surface, precise pressure distribution 
predictions in some chord regions are critical. If the 
shock location is forward of the prediction, or if the 
leading-edge recompression is too steep or extends 
past the suction region, large regions of laminar flow 
could be lost. 

The first calculations performed with the initial 
grid (grid 1) and an early version of USM3D coupled 
to the boundary layer model are compared to the 
experimental data in figures 7 and 8 at Mach numbers 
of 0.793 and 0.820, respectively. Difficulty in accu- 
rately predicting pressures in regions critical to the 
design of a laminar-flow wing is evident at six differ- 
ent chord locations (. x/c locations are specified for the 
52-percent span location) as shown in the table. 


x/c 

Location 

0.00 to 0.04 
0.04 to 0.10 
-0.10 
0.10 to -0.30 
-0.40 
0.40 to 0.50 

Leading-edge favorable pressure gradient 
Leading-edge recompression 
Leading-edge constant pressure 
Midchord favorable pressure gradient 
Midchord shock location 
Midchord recompression after shock 


There may also be an overcompression near the lead- 
ing edge (e.g., fig. 8, 64 percent span). The calcula- 
tions of figures 7 and 8 are clearly inadequate. To 
obtain a better fit between experiment and prediction, 
angle of attack is often changed to match another 
parameter such as lift coefficient. In this work, the 
experimentally measured angle of attack is always 
used. Efforts to improve solution accuracy are docu- 
mented in the following sections. 

6.1. Effect of Experimental Errors 

The maximum experimental errors in Mach num- 
ber and angle of attack were estimated to be ±0.002 
and ±0.25°, respectively. Figure 9 shows the com- 
bined effect of these potential errors by applying the 
positive and negative error extremes to the given flight 
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conditions at M - 0.793 on grid 1 . The two calcula- 
tions do not bracket the experimental data; therefore, 
the disagreement with the experimental data cannot be 
due solely to experimental error. 

A higher Mach number would improve the agree- 
ment between code and experiment. To determine the 
effect of Mach number alone, the M = 0.793 calcula- 
tion was repeated with M = 0.795. The results are 
given in figure 10 where it is seen that this small 
change in Mach number has an almost negligible 
effect on the pressure distribution. 

The M = 0.793, a = 2.598° calculation was 
repeated with a = 2.848° on grid 1 to determine the 
effect of angle-of-attack error alone. The results are 
given in figure 1 1 where it is seen that this change pro- 
vides better agreement in the midchord region but not 
in the leading-edge recompression region (see fig. 11 
inset). On the basis of these and subsequent results on 
a refined grid (appendix A), it is concluded that the 
discrepancy between code and experimental data is 
not due to an angle-of-attack error of 0.25°. 

6.2. Effect of Transition Location 

The previous calculations assume that boundary 
layer transition occurs at the wing leading edge. A 
transition front at 42 to 47 percent chord (the mea- 
sured transition location for the four experimental 
span stations with the suction system active) was spec- 
ified in the input to the code for these stations on 
grid 1. The results are given in figure 12; they show 
that this variation in the exact transition location has a 
negligible effect on the prediction. This issue is revis- 
ited in appendix A on a refined grid. 

6.3. Effect of Boundary Layer 

Inclusion of boundary layer effects can have a 
profound impact on the wing pressure distribution. In 
figure 13, calculations on grid 1 with and without the 
boundary layer model are presented with the experi- 
mental data. As expected, the inviscid solution exhib- 
its a shock farther aft than does the viscous solution. 
The inviscid solution appears to be in slightly better 
agreement with the experimental data than the viscous 
solution. Since all schemes for solving the Euler equa- 
tions contain numerical dissipation that decreases as 
the grid is refined, this agreement is typically an indi- 


cation that the grid is too coarse. The effect of grid 
refinement is investigated in section 6.7. 

6.4. Effect of Nacelle Flow Modeling 

Using a new version of USM3D was necessary to 
overcome a hard-coded limit of 10 boundary layer sta- 
tions in the initial code. This new version also 
included improvements resulting in more physically 
correct nacelle inlet and exit boundary conditions. In 
the original code, the Mach number at the exit of the 
engine is set to the free-stream Mach number, while 
the Mach number at the engine inlet face is adjusted so 
that mass flow is conserved through the engine. The 
newer code uses a more physically correct approach 
(described in ref. 16) in which quantities are set on the 
inlet and exhaust faces of the engine based on the 
directions of characteristic lines at these faces. This 
model allows specification of stagnation conditions in 
the engine and hence simulation of propulsion effects. 
For the cases presented in this work, the conditions 
specified in the engine were those of the free stream; 
therefore, no propulsion effects are simulated. 

Figure 14 shows a comparison of the results 
obtained on grid 1 by using the initial and the 
improved versions of USM3D. By using the older 
nacelle model, the Mach number at the inlet was sub- 
sonic and less than the free-stream Mach number. The 
newer nacelle model yielded a slightly supersonic 
Mach number at the inlet face due to the slight con- 
traction and reexpansion of the diffuser cross section 
upstream of the inlet face. This difference in inlet 
flows has a dramatic impact on the mass flow captured 
by the engine and hence the disturbance of the sur- 
rounding flow field. Both the shock location and the 
leading-edge pressures are significantly improved by 
using the new code, although there is still much room 
for improvement. The improved nacelle boundary 
conditions are used for all remaining calculations on 
grids 2 through 5. 

6.5. Effect of Number of Spanwise Boundary 
Layer Stations 

The number of spanwise boundary layer stations 
needed for an accurate prediction of the experimental 
data was investigated by calculating the wing pressure 
distribution by using 4, 10, 14, 16, and 1 8 stations with 
grid 2 at M = 0.793. As grid 1 was limited to 
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10 span wise boundary layer stations, a new grid 
(grid 2 in table 2) with the capacity for 18 boundary 
layer stations was used in these calculations. The 
results are presented in figure 15. Four stations alone 
are clearly inadequate and result in a pressure distribu- 
tion closer to the inviscid solution. (The inviscid solu- 
tion is not shown.) (Also see fig. 13.) The data 
prediction using 10 stations is nearly as good as the 
prediction using 18 stations (data for 14 and 16 sta- 
tions are indistinguishable from the 18-station calcula- 
tion). We conclude that the number of spanwise 
boundary layer stations is clearly an important factor 
but that the number need not be excessive. 

6.6. Effect of Station Placement 

Figure 16 shows the effect of the placement of the 
spanwise boundary layer stations on the computed 
pressure distribution. In both cases, the number of 
boundary layer stations and the grid resolution are 
fixed. For the calculation denoted as “with experimen- 
tal stations,” the set of boundary layer stations 
included the precise spanwise experimental data sta- 
tions (stations 7, 9, 10, and 1 1 in fig. 5). For the calcu- 
lation denoted as “without experimental stations,” the 
boundary layer data were interpolated to the experi- 
mental data stations. As there is very little difference 
in the computations, we conclude that placement of 
the boundary layer stations is not critical. 

6.7. Effect of Grid Refinement 

The inviscid calculation of figure 13 agreed with 
the experiment better than did the viscous solution, 
and a lack of grid resolution was suspected (see also 
section 6.3). Based on the experimental data, grid 2 
was refined near the leading edge and near the shock, 
yielding grid 3 (table 2 and fig. 2). The impact of this 
refinement is visible in figure 17, which compares 
results on grids 2 and 3. Note the sharper resolution of 
the shock, its better agreement with the experiment, 
and the appearance of a postshock expansion in the 
calculation that agrees well with the experimental 
data. 

The improved agreement with experimental data 
obtained with grid 3 suggested that additional accu- 
racy in the critical leading-edge region could be 
obtained with further grid refinement; however, grid 3 
already required the maximum available computer 


memory. Additional refinement required other regions 
of the grid to be coarsened. Thus, two new grids 
(grids 4 and 5) were generated and are described in 
table 2. Grid 4 further refines the leading-edge region 
of grid 3 while coarsening the fuselage, nacelle, and 
wing trailing edge; grid 5 was generated to assess the 
effects of this coarsening. 

Computed pressures on grid 4 for M = 0.793 are 
shown in figure 18 in comparison with results on 
grid 3. Some improvement of the prediction in the 
leading-edge region is apparent at the 46-percent and 
52-percent span stations. At this level of grid resolu- 
tion, the code is beginning to predict the higher order 
variations of the experimental pressure distribution 
that occur just downstream of the leading-edge peak 
C p . The predicted shock location and pressures down- 
stream of the shock are in excellent agreement with 
the experimental data. The midchord pressure distribu- 
tion prediction is reasonable, but the disagreement 
with experimental data is still not fully explained. 

Figure 19 shows results on grid 5 in comparison 
with those on grid 3 to assess the effects of coarsening 
the fuselage and nacelle. Note that there is no signifi- 
cant loss of accuracy in the regions of primary interest 
(i.e., the region where experimental data are avail- 
able); however, significant loss of accuracy does occur 
near the trailing edge of the airfoil. 

To illustrate the overall effect of the grid refine- 
ments in this study, calculations on grids 2 and 4 are 
presented in figures 20 and 21 for flight Mach num- 
bers of 0.793 and 0.820, respectively. Improvement in 
the agreement of the prediction with experimental data 
is seen in regions where the earlier calculations were 
somewhat inaccurate. 

A comparison of isobar patterns for the two flight 
conditions is also shown in figures 22 and 23. Compu- 
tations are shown on grid 1 by using the original 
nacelle flow model and on grid 4 by using the 
improved nacelle model. Note the relatively parallel 
isobars in the redesigned hybrid laminar-flow wing 
section (located between the two vertical lines) just 
outboard of the nacelle. The strong influence of the 
nacelle pylon on the upper wing surface flow is also 
evident in these figures. 

Appendix B compares the experimental attach- 
ment line location with calculations using grids 1 
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and 4. The results on grid 4 show good agreement with 
the experimental data near the attachment line. 

Based on the sequence of results presented in the 
preceding sections, further refinement of the grid can 
reasonably be expected to produce further improve- 
ments in solution fidelity when sufficient computer 
resources are available. 

7. Concluding Remarks 

Initial calculations did not agree well with experi- 
mental data in six regions of the wing chord. To 
improve calculation accuracy, studies were made of 
discrepancies between computation and experiment 
arising from experimental uncertainties, specification 
of transition location, viscous effects, nacelle flow 
model, number and placement of spanwise boundary 
layer stations, and grid resolution. It was determined 
from these studies that uncertainties in the Mach num- 
ber and angle-of-attack measurements were not signif- 
icant contributors to the initial disagreement between 
calculations and experimental data. It was also deter- 
mined that 10 boundary layer stations were sufficient 
to capture the pressure distribution. 

Availability of flight data was invaluable in guid- 
ing the calculations. The first refinement of the in vis- 
cid grid (from grid 2 to grid 3) resulted in very little 
change in the pressures in the leading-edge region, and 
it is unlikely that further refinement would have been 
done solely on the basis of the computed results. The 
experimental data, however, indicated subtle varia- 
tions in the leading-edge surface pressures that the 
calculation on the refined grid (grid 3) did not capture. 


On this basis, the second refinement (from grid 3 to 
grid 4) was performed. 

Improvement of the nacelle flow model had the 
largest single impact on the accuracy of the pressure 
distributions. Of particular importance is the 
10-percent chord rearward movement of the shock due 
to the change in nacelle flow modeling. The upper sur- 
face rooftop Mach number increased by approxi- 
mately 0.05, and the pressures near the leading edge 
were greatly improved. 

Refinement of the inviscid grid also had a signifi- 
cant impact on accuracy. Two refinements of the 
inviscid grid concentrating primarily on the leading- 
edge region resulted in significantly improved calcula- 
tions of the leading-edge pressures, the shock location 
and definition, and the pressures downstream of the 
shock. Reduction of the grid spacing at the leading 
edge by approximately 70 percent to 0.05 percent of 
the mean aerodynamic chord (MAC) was required to 
achieve good agreement between the calculations and 
experiment in five of six chord regions where the 
initial calculations were in poor agreement. The 
improved grids, however, were unable to capture the 
small local pressure gradients occurring in the leading- 
edge region of the outboard data stations. Further grid 
refinement is required to obtain these details; 
however, the limits of current computer resources 
were used in this study. A solution on a grid of 
1 250233 cells (grid 4) required 256 million words of 
memory and approximately six hours of processing 
time on the Cray C90 computer. 

NASA Langley Research Center 
Hampton, VA 23681-2199 
April 1, 1998 
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Appendix A 

Effect of Angle-of- Attack Errors, 
Transition Location, and Viscous 
Effects on a Refined Grid 

Several calculations were repeated for the 
M = 0.793 case on grid 4 to confirm conclusions 
drawn from results obtained on grid 1 . 

The first case represents an increment in angle of 
attack of 0.25° to determine the effect of a maximum 
angle-of-attack error. The results are given in 
figure 24 and show that the increased angle-of-attack 
prediction is worse at the leading edge and that the 
shock is predicted to be aft of its true location. This 
result confirms the conclusion drawn in section 6.1 
that discrepancies between experimental data and 
computations were not due to an experimental error in 
angle-of-attack measurement. 


Figure 25 shows the results obtained by specifying 
a transition front at 42 to 47 percent chord in the mid- 
span region. Although the corresponding calculation 
on grid 1 (fig. 12) was virtually identical to the fully 
turbulent result, the transitional case on grid 4 shows 
less agreement with the experiment than with the fully 
turbulent case. Because the experimental data were 
taken with the flow fully turbulent, it is expected that 
the fully turbulent calculation will provide the best 
prediction. 

In figure 26, calculations with and without the 
boundary layer model are presented for grid 4. The 
same behavior noted in section 6.3 (fig. 13) is present 
in that the inviscid solution exhibits a shock well aft of 
the shock location of the viscous solution. However, in 
figure 26, the viscous result is in much better agree- 
ment with the experiment than is the inviscid result. 
The reduction of numerical dissipation due to grid 
refinement has moved the shock aft in both the invis- 
cid and viscous solutions shown in figure 26. 
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Appendix B 

Comparison of Experimental and 
Calculated Attachment Line Location 
for Two Grids 

Detailed study of wing boundary layer flow, 
including transition prediction and stability analysis, 
often requires nearly exact knowledge of the attach- 
ment line location. It is therefore of interest to com- 
pare experimental and analytical results for the 


attachment line and also to see whether and how the 
grid size used might affect the location of the attach- 
ment line. Results in figures 27 and 28 for the four 
span locations of this study show that grid 4 does a 
better job of predicting the upper surface pressures 
near the highlight. The experimental attachment line 
location for the present conditions is near the highlight 
and agrees well with the grid 4 prediction. The pre- 
dicted attachment line location is within 0.4 in. for the 
two grid sizes used in this study (table 2), but it is 
noted that the predicted attachment line may fall on 
opposite sides of the highlight, depending on the grid. 
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Table 1 . Summary of Flight Conditions and Geometry for Experimental Data 


Parameter 

Case 1 

Case 2 

Mach number 

0.793 

0.820 

Angle of attack, deg 

2.598 

2.377 

Altitude, ft 

34417 

37 176 

Reynolds number based on MAC 

31.9 x 10 6 

30.0 x 10 6 

Lift coefficient 

0.447 

0.449 

Sonic pressure coefficient 

-0.455 

- 0.379 

Span stations, percent 

46, 52, 57, 64 

Semispan, in 


750 

MAC, in 


200 


Table 2. Summary of Grids Used in Calculations 



Global statistics 

Average spacing in test region, 
percent MAC, at nominal xic — 

Description 

Cells 

Nodes 

(total) 

Boundary 

nodes 

0 

0.15 

0.5 


HOSH 

891 987 

164 876 

25 388 



MEM 

0.836 

Initial grid, 10 IBL stations 

2 

895 759 

165 647 

25 618 

1 


mm 

0.836 

Expanded to 1 8 IBL stations 

3 

1 341 188 

244 211 

30079 

0.134 

0.251 

1.338 

0.836 

Refined leading-edge and 
midchord region 

4 

1 250233 

229 052 

31 724 

0.050 

0.084 

1.338 

1.672 

More leading-edge refinement, 
coarsened fuselage and 
nacelle 

5 

598 243 

110510 

16 761 

0.134 

0.251 

1.338 

1.672 

Coarsened fuselage and nacelle 
without leading-edge 
refinement 
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(b) Oblique view. 

Figure 1. Leading-edge test article and pressure belt locations on wing surface. 
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(c) Grid 4. 


Figure 2. Test surface grid refinement. 
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(a) Aircraft surface. 



(b) Wing and nacelle. 

Figure 3. Surface triangulation of aircraft geometry (grid 4). 
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Figure 4. Boundary layer station locations for grid 1. 




Figure 6. Pressure distributions typical of laminar-flow wing. 
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Figure 11. Effect of potential experimental angle-of-attack error on surface pressure coefficient distribution on grid 1 
(M = 0.793; a = 2.598°; R = 31.9 x 10 6 ; 10 EBL stations). 
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Figure 14. Comparison of surface pressure coefficient with experiment for results from two versions of nacelle inflow/outflow 
boundary conditions on grid 1 ( M = 0.793; a = 2.598°; R = 31.9 x 10 6 ; 10 IBL stations). 
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distribution measurements. 
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distribution measurements. 
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